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ABSTRACT 

Recent developments in numerical relativity have made it possible to follow 
reliably the coalescence of two black holes from near the innermost stable circular 
orbit to final ringdown. This opens up a wide variety of exciting astrophysical 
applications of these simulations. Chief among these is the net kick received 
when two unequal mass or spinning black holes merge. The magnitude of this 
kick has bearing on the production and growth of supermassive black holes during 
the epoch of structure formation, and on the retention of black holes in stellar 
clusters. Here we report the first accurate numerical calculation of this kick, for 
two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on 
analytic considerations to give a significant fraction of the maximum possible 
recoil. Our estimated kick is 105 km s -1 with an error of less than 10%. This 
is intermediate between the estimates from two recent post-Newtonian analyses 
and suggests that at redshifts z > 10, halos with masses < 10 9 M© will have 
difficulty retaining coalesced black holes afer major mergers. 

Subject headings: gravitational waves — relativity 
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1. Introduction 

When two black holes merge, the gravitational waves they produce will carry away net 
linear momentum, barring conditions of special symmetry (e.g., two equal-mass nonspinning 
black holes receive no kick). The magnitude of this recoil is important in a variety of 
astrophysical contexts, including the cosmological evolution of supermassive black holes 
(Merritt et al. 2004; Boylan-Kolchin, Ma, k Quataert 2004; Haiman 2004; Madau k 
Quataert 2004; Yoo k Miralda-Escude 2004; Volonteri k Perna 2005; Libeskind et al. 2005; 
Micic, Abel, k Sigurdsson 2005) and the growth and retention of intermediate-mass black 
holes in dense stellar clusters (Taniguchi et al. 2000; Miller k Hamilton 2002a, b; Mouri 
k Taniguchi 2002a, b; Miller k Colbert 2004; Gultekin, Miller, k Hamilton 2004, 2006; 
O’Leary et al. 2006). There is therefore a long history of analytical estimates of this recoil 
(Peres 1962; Bekenstein 1973; Fitchett 1983; Fitchett k Detweiler 1984; Redmount k Rees 
1989; Wiseman 1992; Favata, Hughes, k Holz 2004; Blanchet, Qusailah, k Will 2005; 
Damour k Gopakumar 2006). However, it has been shown that almost all of the recoil 
occurs in the strong gravity regime, inside the innermost stable circular orbit (ISCO). This 
is precisely where analytical treatments are least reliable. An accurate estimate of the recoil 
kick therefore requires full numerical simulations of the final phase of the inspiral, merger, 
and ringdown of the coalescence of two black holes. 

Until recently, numerical simulations were not stable and accurate enough for such 
estimates. This situation has changed dramatically in the past year, with several groups 
developing codes that have allowed the evolution of binary black hole spacetimes from 
close to the ISCO through merger and ringdown. These codes utilize techniques which 
allow the black holes, with their inherent singularities, to move successfully through the 
computational domain. One approach (Pretorius 2005; Pretorius 2005) is to excise the 
singular region in the (physically inaccessible) interior of the black hole. In contrast, our 
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method allows the singular region to be approximately represented in the computational 
domain (Bakeret al. 2005, 2006; Campanelliet al. 2005b). These new techniques have led 
to dramatically more effective numerical simulations of binary black hole systems, recently 
allowing us to determinate accurate waveforms for the final orbits and merger of equal 
mass systems (Bakeret al. 2006), which we generalize here for nonequal masses. While 
some preliminary numerical calculations of recoil from mergers of non-spinning black holes 
have been reported (Campanelli 2005a; Herrman, Shoemaker, & Laguna 2006), the initial 
separations have been too small and the resolutions too coarse for reliable and precise 
numbers. 

Here we report the first precise fully numerical estimates of the kick received from the 
merger of two nonspinning black holes. We choose a 1.5:1 mass ratio because it is close to 
the analytically estimated optimal mass ratio for maximum kick, but is also close enough to 
equal mass that resolution issues are not serious impediments to the numerical evolution. In 
§ 2 we describe our numerical method and results, including convergence tests. We discuss 
the astrophysical implications of these results in § 3. 

2. Numerical Simulations 

Reliable simulations of binary black hole mergers requires the specification of 
constraint-satisfying initial field data followed by stable evolution of the Einstein equations. 
The evolution variables hold information about gravitational fields in the form of tensor 
fields representing the curvature of a vacuum spacetime. Einstein’s equations, together with 
a crucial specification of gauge (i.e., coordinate) conditions, then govern the evolution of 
these fields in time. 

We utilize the “puncture” approach (Brandt & Briigmann 1997) to specify constraint- 
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satisfying initial field configurations for black holes on an approximately circular inspiral 
trajectory within an orbit or two of merger. Specifically, we consider two cases of initial 
configurations, with the black holes at initial separations of dj n ; t = 4.1 M 0 and d in it = 6.2 M 0 , 
where Mo is the total initial gravitational mass of the system, and we use units where 
G = c = 1 so that all quantities can be represented in terms of their mass-scaling. 


Hahndol (Imbiriba et al. 2004), using a second-order iterative Crank-Nicholson scheme, 
on an adaptive mesh refinement structure implemented via PARAMESH (MacNeice et al. 
2000). The punctures are allowed to move freely through the grid and are evolved according 
to a variant of theJBSSN formulation of Einstein’s equations and gauge conditions as 
described in BaWret ail. (2006). 


We interpret our numerical results by studying the gravitational radiation generated by 
the merger which carries gauge-invariant information away from the strong-field region. In 
our simulations, radiation is represented by a component of the space-time Weyl curvature 
tensor, ?/> 4 . In terms of ipi, the net momentum radiated by the system can be expressed as 
follows: 


dPi 


= lim 


<m- 

r 


^4 


( 1 ) 


dt r-ioo ^ 47T 

We compute ?/> 4 numerically, and extract it from the simulation data on a sphere of radius 
r = 40M, which we have found is sufficiently large enough to give accurate results. 


Our main numerical results are given in Figure 1, which shows the net speed of the 
center of mass as a function of time for several runs. To compare runs at different initial 
coordinate separations we have shifted the time axis so that the kick speed peaks at t = 0. 
Note that for the d in ; t = 4.1Mo case we show the results of two runs with different numerical 
resolutions, which give kick histories and final kicks that are almost indistinguishable, 
differing by less than 1%. This provides an indication of the accuracy of the simulations. 
Note also that the run with separation d- in i t = 6.2 M 0 gives a final kick just 10% less than the 
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closer runs, suggesting that the sensitivity of the result to the choice of initial configuration 
has only a small effect, which we expect to be less important for the case of larger initial 
separation. 

The inspiral from infinity to our starting point would in reality have imparted a kick 
to the center of mass at the beginning of our runs. To explore the effects of this in a simple 
way, we added an initial momentum vector to our larger separation run, with a magnitude 
of a few km s -1 as suggested by post-Newtonian analyses (e.g., Blanchet, Qusailah, & Will 
2005) and with a direction chosen to smooth out the evolution of the kick. The result is 
shown in the solid line in Figure 1, and ends up with almost exactly the same final kick. 
This will need to be confirmed by runs with larger initial separation, but it does suggest 
that, as expected analytically, the net kick grows smoothly and monotonically until the 
merger. 

Recently, two groups have made refined analytic estimates of the kick from merger of 
nonspinning black holes, both with precise answers but differing in magnitude by a factor 
of three. Blanchet, Qusailah, & Will (2005) predict a speed of 155±25 km s -1 for a 1.5:1 
mass ratio, whereas Damour & Gopakumar (2006) predict 48 km s -1 . Our result is between 
these estimates but inconsistent with either one. An initial numerical estimate for 1.04:1 
and 1.18:1 mass ratios was made by Herrman, Shoemaker, & Laguna (2006) (albeit at much 
lower resolution and starting much closer than our simulations), and they find kick speeds 
of 9 km s' 1 and 33 km s -1 , respectively. One can compare results at different mass with 
reference to the Fitchett (1983) fitting formula kick oc q 2 (q — 1) /(I + q) 5 = (/ i/m) 2 dm/m , 
where q = mi/ m 2 > 1 is the mass ratio, \x/m — mim 2 /(mi + m 2 ) 2 is the symmetric mass 
ratio and dm/m = (mi — m 2 ) /{mi + m 2 ) is the fractional mass difference. Although 
originally derived in the context of a leading order post-Newtonian approximation, this 
formula is suggested to be a good approximation by perturbative Schwarzschild calculations 
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(Fitchett & Detweiler 1984), and it also closely agrees with more recent post-Newtonian 
results (Blanchet, Qusailah, & Will 2005; Damour & Gopakumar 2006). By Fitchett’s 
formula, the Fterrman, Shoemaker, & Laguna (2006) results would imply ~ 80 km s _1 for 
a mass ratio of 1.5:1, closer to our value than the analytic estimates but still significantly 
below ours. However, the convergence of our results in the numerical and physical 
domains suggests that our kick is correct to < 10% and can therefore be used to test 
strong-field extrapolations of post-Newtonian analyses as well as to investigate astrophysical 
consequences. 


3. Discussion and Conclusions 

Our estimate of 105 ± 10 km s _1 for the kick received in a 1.5:1 merger of nonspinning 
black holes has important implications for the assembly of supermassive black holes in the 
early universe (see also Merritt et al. 2004). This is because as dark matter halos merge in 
the process of hierarchical structure assembly, their central black holes are also presumed to 
merge if the halo mass ratio is not too extreme (otherwise dynamical friction on the halos is 
likely to be inefficient; see Taffoni et al. 2003). If black hole masses are linked to the mass 
of their host halos, it is therefore expected that in the early universe, black hole mergers are 
likely to involve comparable-mass objects. If the resulting kick exceeds the escape speed of 
the merged dark matter halo, the halo is left without a black hole. This could, therefore, 
have a significant impact on the number of black hole mergers in the early universe. In 
addition, if mergers between comparable-mass halos are common at redshifts z > 10, these 
ejections might reduce substantially the fraction of halos that host black holes. 

A full appraisal of the consequences of recoil will require detailed numerical studies of 
the effects of black hole spin and mass ratio, which will be our subject in future papers. 
We can make initial guesses by focusing on nonspinning black holes and by adopting as 
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before the Fitchett (1983) formula for the dependence of kick speed on mass ratio. Setting 
u kick (g = 1-5) = 105 km s _1 fixes the curve. 

To estimate the escape speed from a dark matter halo, we follow the treatment of Merritt 
et al. (2004). The escape speed from a halo of virial mass M v , r is U e 2 c = 2cg(c)M yiT /R V i T , 
where as before we set G = 1. Here c is the concentration parameter for the halo, 
g(c) = [ln(l + c) - c/(l + c)}~\ and we follow Navarro, Frenk, & White (1997) in defining 
Ry ir as the radius inside which the average density is 200 times critical (and thus ignoring 
the very weak dependence on cosmology): Ry ,- r = [M v ; r /( 200 x 47rp cr j t /3)] 1//3 , where 
p CT it = 31T 2 /(87r) « 10~ 29 g cm~ 3 (l + z) z . Bullock et al. (2001) find that the concentration 
parameter scales as c = 9 [M vir /(2 x 10 13 M©)] -0 ' 13 (1 + z)~ x . The resulting escape speed is 

F eS c ~ 220 km s _1 [^(c)] 1/2 (M vir /10 10 M©)°‘ 27 . (2) 

In Figure 2 we plot the minimum halo mass as a function of redshift such that 
V esc > Ukick for the listed mass ratios of 1.5:1, 3:1, 5:1, and 10:1. Note that the low power 
of My ir in the escape speed means that the minimum halo mass depends sensitively on the 
kick speed. For example, if the kick speed at a 1.5:1 mass ratio were 200 km s -1 instead 
of 105 km s _1 , the threshold halo mass at z = 10 for this mass ratio would jump from 
4 x 10 8 Mq to 3 x 10 9 M q . This would in turn reduce the number density of halos of at 
least the required mass by a factor ~ 20 (see Mo & White 2002 for a pedagogical discussion 
of how to estimate halo number densities). This points out the importance of estimating 
kick speeds accurately. 

As discussed by Merritt et al. (2004), kicks are also important in the current-day 
universe for low-mass concentrations of stars. Their Figure 2 is a useful summary of escape 
speeds from the centers of galaxies and globular clusters. From their figure, we see that 
comparable-mass mergers with kick speeds ~ 100 km s -1 will cause ejection from globulars 
and some dwarf galaxies, but that more massive galaxies will retain the remnant. Using 
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the Fitchett (1983) scaling with mass ratio, we find that black holes that are > 10 times 
more massive than their companions are not ejected from globular clusters with escape 
speeds of a few tens of km s -1 . Coincidentally, this is approximately the same mass ratio 
required to protect a massive black hole from cluster ejection from three-body interactions 
(Giiltekin, Miller, k Hamilton 2004). Therefore, if an intermediate-mass black hole of 
> 100 — 200 Mq is formed in a cluster, it can stay and potentially grow through future 
interactions. Stellar-mass binary black holes (M < 50 M Q ) will be ejected from such clusters 
prior to merger by Newtonian three-body interactions (Kulkarni, Hut, k McMillan 1993; 
Sigurdsson k Hernquist 1993), hence recoil is not so important for low-mass black holes in 
this context. However, mergers of black hole binaries in low-density galactic disks could 
produce a population of high-speed coalesced black holes. 

In conclusion, we have presented the first reliable and precise fully numerical estimate 
of the gravitational recoil produced by the merger of two unequal mass nonspinning black 
holes. Our value of 105 km s -1 for a 1.5:1 mass ratio is accurate to within < 10%, as 
determined by our runs at different resolutions and different initial separations. This speed 
is intermediate between the recent analytic estimates of Blanchet, Qusailah, k Will (2005), 
who suggest 155±25 km s -1 for this mass ratio, and Damour k Gopakumar (2006), whose 
formulae would imply 48 km s _1 . Our results are thus an important first step in the 
evaluation of the astrophysical consequences of gravitational radiation recoil in dense stellar 
clusters and the early universe. 
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Programs Office of the National Research Council and Oak Ridge Associated Universities. 
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Fig. 1. — The magnitude of the radiated momentum, as a function of time, from three 
different simulations. For an initial separation between the black hole punctures of coordinate 
distance d init = 4.1 M 0 , the final values of the momenta from two different resolutions are 
shown to agree to within 1%. For an initial separation of coordinate distance d init = 6.2M 0 , 
the final value of the momentum is shown to agree with that of the closer initial separation 
to within 10%. Also shown, the dj n ; t = 6.2M 0 momentum is recomputed assuming a nonzero 
initial momentum vector P in j t , estimated so as to give a more monotonically increasing curve 
and presumably thereby to approximate the momentum radiated up to this point by black 
holes of initially infinite coordinate separation. The final kick is nearly the same for all 
curves, strongly suggesting that the result is correct to better than 10%. 
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Fig. 2. — Minimum mass of a dark matter halo at a given redshift required to retain the 
product of the merger of two nonspinning black holes with a mass ratio indicated on the curve. 
Details of the computation are in the text; note that we use the Fitchett (1983) analytical 
estimate of the mass ratio dependence. This figure indicates that early halos might lose 
merger remnants because of the kick, but in the current universe only the smallest dwarf 
galaxies or globular clusters could have black holes ejected because of gravitational radiation 
recoil. 



